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The Lichnerowicz and Israel theorems are extended to higher order theories of 
gravity. In particular it is shown that Schwarzschild is the unique spherically sym- 
metric, static, asymptotically flat, black-hole solution, provided the spatial curvature 
is less than the quantum gravity scale outside the horizon. It is then shown that in 
<^ • the presence of matter (satisfying certain positivity requirements), the only static 

and asymptotically flat solutions of General Relativity that are also solutions of 
k*" | higher order gravity are the vacuum solutions. 
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O ' I. INTRODUCTION 

cr 

Whilst General Relativity (GR) remains the most successful classical theory we have, 
one may question whether it is possible for there to be corrections to it. Indeed we expect 
GR to break down as we approach the Planck scale where the quantum theory of gravity 
(whatever it is) will become dominant. In this sense GR is not a full physical theory, only 
an excellent approximation to some (presumably) complicated underlying theory, with the 
approximation getting better and better as the scales we consider are further and further 
from the Planck scale 1 . 

Over the past century there have been many important theorems proved about the so- 
lutions and structure of GR and an important question is whether these theorems are also 
valid in theories that deviate from GR at some scale. In this paper we will look in particular 
at two theorems, firstly the Lichnerowicz theorem [2], which tells us that the only static, 
asymptotically flat, geodesically complete, vacuum, solution to Einstein's equations is flat 
space-time. We then consider the Israel theorem which demonstrates that the only static, 
asymptotically flat, vacuum space-time, which contains past and future event horizons (that 
intersect on a surface that is topologically S 2 ) is given by the Schwarzschild metric 2 . This 
'no-hair' theorem is a striking result in classical GR and extensions of it to more general 
gravity theories, in particular as we approach the scales on which GR is expected to be 
violated, would provide us with insight into the transition from GR to full quantum gravity. 

The action for GR is the well-known Einstein-Hilbert action which is linear in the Ricci 
scalar R (here for simplicity we neglect the cosmological constant). A natural extension of 
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1 It is important to note that here the Planck scale can include both very small (UV) scales and very large 
(IR) scales (see for example [lj]). 

2 The Israel theorem also shows that the Reissner-Nordstrom metric must be the solution if the vacuum is 

replaced by an electromagnetic field. 
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this is to include terms of higher powers of curvature invariants to find, 

<S Grav . = J d 4 xv^ [^R - a hR abcd R abcd - ahR ab R ab + (3hR 2 + O (i? 3 )) , (1.1) 

where we have written only terms that are only 2 nd order in the curvature. Through-out 
we will use Latin indices to label space-time components a, b, ■ ■ ■ = 0, . . . , 3, our signature 
is (— , +, +, +) and our conventions for curvature are R abc d = d b T d ac + . . . and R ab = R c acb - 
In Eq. (11.11) the coefficients ao, a, (3 and 7 are dimensionless number and k 2 = 32ttG, 
so GR is recovered simply by setting the coefficients of the higher order terms to zero 
(«o = ot = (3 = 0) and taking 7 = 2. Also note that the presence of the h in Eq. (11.11) is only 
due to dimensions i.e. this remains a classical theory. If these correction terms are motivated 
from some quantum gravity theory, then we may expect the coefficients ao, a and (3 to be 
order unity, however in general they can take any value. 

In fact, in 4-dimensions, the Gauss-Bonnet combination of second order curvature invari- 
ants integrates to a topological invariant, i.e. 

R abcd Rabcd ~ 4:R ab R ab + R 2 = total divergence , (1.2) 

which we can use to eliminate the R abcd R abca - term from Eq. (11.11) . Thus, up to second order 
in the curvature, the most general correction to GR comes from the action, 

<S Grav . = J d A x^g~ (^R - ahR ab R ab + f3HR 2 ) . (1.3) 

This is the theory that we will work with and in the following we refer to it as 4 th order 
gravity, since the equations of motion involve, at most, 4 th derivatives of the metric. Histor- 
ically, theories of this type, with 7 = were considered for many reasons, not least because 
of their improved behaviour under renormalisation, however it was shown that without a 
term linear in R, the theory (for a = 0) does not couple to matter correctly 0]. The 
theory, with general 7 was first considered in [5| and later in j6[, where a perturbative anal- 
ysis around the Schwarzschild solution was performed, whilst the initial value formulation 
of the theory was given in Q. The consequence for black holes in higher dimensions has 
been consider in [gj and the cosmological implications of 4 th order gravity were investigated 
in 0-12], with the extension to include (tensor) perturbations appearing only recently (l3| . 
The special case in which the correction terms are conformal (when a = 3/3, which gives 
the Bach- Einstein equations) has received considerable attention (see for example [13]) and 
has recently appeared in the context of non-commutative geometry [TBI], with much work 
going into restricting the par ameters of this form of the theory via cosmological and astro- 
physical observations [16|-|20|. More gen erally, corrections to GR of this type are expected 
from String Theory (see for example |2l[|). renormalization group techniques j22| and Loop 
Quantum Gravity [231 ] . 

In addition there has been a great deal of phenomenological work on the cosmolog- 
ical consequences of higher order gravity theories, usually in the form of / (R) theories 



see for example [24H26I]). Despite this it is important to note that theories defined by 



Eq. (11.31) cannot represent true physical theories for arbitrary coefficients due to the pres- 



ence of ghosts [6|, |27|, [28| . Not only are ghost a problem for the quantisation of the theory, 



they also lead to a break-down of causality at the classical level. Note however that with 
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certain restrictions on the coefficients (e.g. a = and < 0), all the difficulties associated 
with these ghosts can be eliminated. Here we will take the (conservative) view that, just as 
for GR, 4 th order gravity should be considered only as an approximation to the underlying, 
full theory (which is presumable free of ghosts etc.). 

In this paper we shall demonstrate three key results, the first two: an extension of the 
Lichnerowicz (Section (IIIII) ) and Israel (Section ( II VI) ) theorems to 4 th order gravity (for 
space-times satisfying some restriction on the magnitude of the spatial curvature) suffice to 
demonstrate that Schwarzschild is the unique, spherically symmetric, static, asymptotically 
flat, vacuum solution to the theory 3 . The final result (Section (jyj) considers non-vacuum 
solutions to the theory and demonstrates that for static, asymptotically flat space-times 
with matter satisfying certain positivity requirements, it is only the vacuum solutions of GR 
and 4 th order gravity that agree; all other solutions of GR fail to be solutions of 4 th order 
gravity. 



Considering the vacuum, the trace part gives the Klein-Gordon equation for the Ricci scalar, 



It is well known that there are no static, asymptotically constant solutions of the Klein- 
Gordon equation for a scalar field (see for example [29']), however since the method used 
to prove this result will be the one that is later extended to the full equation of motion, 
Eq. (12. ip . we give it explicitly here. 

Theorem II. 1. If the space-time is static (C t g a b — 0j and D a R — > sufficiently fast at 
infinity, then (□ — m?) R = 0, for m ^ G 1R implies R = 0. // the space-time is static and 
R — > sufficiently fast at infinity the equation UR = implies R = 0. 

Proof. We decompose the metric into components parallel and perpendicular to the constant 
time space-like hyper-surfaces S via g a b = h a b — \t a t bl where t a is the time-like Killing field, 
which has norm t a t a = —A and h ab is spatial metric on S. Note the following useful identities 



II. STATIC SOLUTIONS OF THE EQUATIONS OF MOTION 




^T ab . (2.1) 




(2.2) 



C t R = t b S7 b R = , 



(2.3) 



1 



(t b V a X - t a v b 



A) 



(2.4) 



6 2A 



3 This result holds for black-holes in which the spatial curvature outside the horizon is small relative to 
the scale set between jk~ 2 and ha. If these corrections come from some quantum gravity theory, this 
restriction is essentially that the spatial curvature be small compared to the quantum gravity scale. 
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and 

t a V a X = . (2.5) 
Using these one can readily show that, 

g ab V a V b R - m 2 R = D a D a R + — (D a X) (D a R) - m 2 R = , (2.6) 

2X 

where D a is the spatial covariant derivative compatible with h ab . Multiplying this equation 
by \ l / 2 R and integrating it over the entire spatial slice, one finds 



/ Vhd 2 x X 1/2 RD a D a R + -\-^ 2 R (D a X) (D a R) - m 2 X 1/2 R 2 
Js L 2 

Integrating the first term by parts this becomes, 

J Vhd 3 x [D a {X 1/2 RD a R) - X 1 ' 2 (D a R) (D a R) - m'X^R 2 } = . 



. (2.7) 



(2.8) 



If the spatial slice is asymptotically constant i.e. D a R — > sufficiently fast at infinity, then 
the boundary term vanishes. The second and third terms in the integrand of Eq. (12. 8 p are 
manifestly negative definite and hence each term must vanish at every spatial point. If m ^ 
this implies R = 0, for the case of m = this gives D a R = and hence we additionally 
require that R — > at infinity in order to find R = 0. □ 

Note that here we require only that D a R — > (or R — > 0) sufficiently fast at infinity, 
which is satisfied by asymptotically constant (or flat) space-times, but is in general a weaker 
condition. In the following we will consider asymptotically constant (or flat) space-times 
since they are of more interest, however all the results hold also for the weaker condition 
given above. 

The basic method of this proof is the following. Decompose the metric onto (and perpen- 
dicular to) constant time space-like hyper-surfaces and use this to write the 4-dimensional 
covariant derivatives as spatial, 3-dimensional covariant derivatives plus correction terms. 
Multiply the resulting equation by a suitable factor and integrate over the space-like hyper- 
surface. Prove that, up to boundary terms, the integrand has a definite sign and hence 
vanishes at each spatial point. This is precisely the same method that we will use to prove 
that, with some restrictions on the magnitude of the spatial curvature on S, there are no 
static, asymptotically constant solutions to the trace-free part of Eq. (12. ip . which, using 
theorem flll.ip can be written as, 

ha ^DR a b — ^ ab 2^ ab ^ Cd ^' cd ~ 2-R cd -Racfe(i^ = , (2.9) 

for a — 3(3 > 0. Since Eq. (I2.9P is rather more complicated than Eq. (I2.2p . we proceed by 
first proving a series of lemmas 4 . 

Lemma II. 2. If the space-time is static and asymptotically constant, then (□ — m 2 ) R ab = 0, 
for m 7^ G E implies R ab = 0. // the space-time is static and asymptotically fiat, then 



A similar line of argument has used in j^sj . 



nR ab = 0, implies R ab = 0. 

Proof. We begin by noting several identities, which are proved in Appendix IVIl 

t a V a t b V b R cd = 2 [{R a{c V d) t b ) (V b t a ) + R ab (V c t b ) (V d t a )] (2.10) 

(V d t c ) (V a t d ) = ^t c t a (V d A) (V d X) + -L (V C A) (V a A) , (2.11) 

^h ab (V a A) t c t d V b R cd = i (V ffl A) (Vfi?) , (2.12) 

and 

jt c t d h ab V a h e b V e R cd = D a D^R , (2.13) 

where ^R is the 3-dimensional Ricci scalar formed from the spatial metric h ab . Now, 
decomposing g ab = h ab — jt a t b one finds, 

(□ - m 2 ) R ab = h ab V a h e b V e R cd - jt a t b V a V b R cd - m 2 R cd = . (2.14) 

Noting the staticity of R ab , 

C t R ab = t a V a R cd + R a <N c t a + R ca V d t a = , (2.15) 
and using Eq. ( 12.10D and Eq. ( 12. lip , one can show (see Appendix IVIl) that Eq. ( 12. 14ft becomes, 

h ab V a h e b V e R cd + ^h ab (V a A) (V feJ R cd ) + ^ [t c t a (V 6 A) (V 6 A) - A (V C A) (V a A)] i? ad 

[t a V,A - t d V a A] [t fe V c A - t c V d A] R ba + ^ [t d t a (V 6 A) (V fe A) - A (V d A) (V a A)] R ac 
-m 2 R cd = . (2.16) 

We now project this equation onto the spatial slice, S with h c e h d and define the projection 
of R nb onto S as R nb = h c „hfR rd to find, 



D a D a R ef + i- (D a X) {D a R ef ) - (D e X) (D a X) R af - (L> /A ) (D e A) t a t fe i? afe 
' ( J D a A)( J D / A)i? ea -m 2 i? e/ = . (2.17) 
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Noting that the extrinsic curvature of the spatial slice vanishes (as can be shown by direct 
calculation) the Codacci and Gauss equations reduce to, 



h b a t c R bc = , (2.18) 
1 t a t b R ab = ^R, (2.19) 



A 



where we have use Theorem (IH.lj) to set R = 0. Using these relations in Eq. ( I2.17P we find, 



D a D a R cd + -L (D a X) (D a R cd ) - ^ (D C X) (D a X) R ad - ^ (D d X) (D C X) {3) R 
1 



{ . 2 (D a X) (D d X) R ca - m 2 R cd = . 



(2.20) 



Multiplying this by X 1 ^ 2 R cd and integrating over the 3-dimensional spatial slice gives, 
J y/hd 3 x D a (X^ 2 R cd D a R cd ) - X 1 ' 2 (D a R cd ) (D a R cd ) - ^X~ 3/2 R cd R ad (D C X) (D a X) 

. (2.21) 



1 



__ A -3/2 (3) RR cd (DdA) p cA) _ x 1 / 2 m 2 R cd R cd 



We now project Eq. (12. 16[) perpendicular to the space-like hyper-surface, with t b t c and 
again use the Codacci and Gauss equations (Eq. (I2.18P and Eq. (I2.19P respectively) to find, 

7 6\\(3) 



tH d h ab V a h e b V e R cd + ^t c t d h ab (V Q A) V b R cd - ^ (V 6 A) (V b A) w R 

~^x (yaA) ( V ' A ) Rab ~ Xm2 {3)R = - 



(2.22) 



Multiplying this by A 1 ^ 2 ^ R, integrating over the 3-dimensional hyper-surface and using 
the identities given in Eq. ( I2.12P and Eq. ( 12 .131) . gives, 



Vhd 3 x D a (A 1 / 2 ^ RD^R) - X 1/2 (D a ^ R) (D a ® R) 
l -X- z ' 2 (D a X) (D a X) (^R) 2 - V 3 / 2 (D a X) (D b X) R ab ®R - X 1/2 m 2 (®R)' 



= . 
(2.23) 



Finally, adding Eq. fl2T2T|) and Eq. fl2~23|) one finds, 



Vhd 3 x 



D a (\V2 R cd DaRcd + A V2 0) RDa (3)^ _ yl/2 (pa R cd^ (/^^ 



1 



_ A l/2 ^ D a (3) R j (3)^ _ _ A -3/2 ^cd^ + (3) jyjd^ (^/^ + (3) RD d X) 



X 1/2 m 2 (R cd R cd +(^R) [ 







(2.24) 



This is sufficient to prove the desired result, however there is a simplification that can be 
made by noting that the contracted Bianchi identity V a (R a b — \dabR) = implies that 



R dc D d X + (3) RD C X = -2XD a R ac . 



(2.25) 
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Substituting this into Eq. (I2.24p explicitly removes the derivatives of A from the integrand 
and ensures that A only appears with positive powers. This will be crucial when we consider 
space-times with (null) horizons. 

In any case, asymptotic constancy i.e. D a R cc i — > at infinity, ensures that the boundary 
terms in Eq. C I2.24I) vanish. All the remaining terms are negative definite and hence each 
term must independently vanish at every spatial point. For m 7^ this implies = 
and R c d = 0, which together with theorem ( 111.11) and Eq. (12. 181) implies R a b = 0. For the 

(3) — 

m = case we find that D a R = and D a R c d = and hence additional require asymptotic 
flatness in order to find R a b = 0. 

□ 

This result is easily extended to include the third term in Eq. ( 12. 9 p by using the fact that, 
by Theorem ( 111.11) we can set R = 0. 

Lemma II. 3. If the space-time static and asymptotically constant, then the pair of equations 
R = and (□ — m\) R ab + m 2 g ab R cd R c d = 0, for mi ^ G R and m 2 G R imply R ab = 0, 
whilst for rrii = 0, additionally requiring the space-time to be asymptotically flat implies 

Rab = 0. 

Proof. From the equation (□ — ml) R ab + m\g ab R cd Rai = it follows that, 

[\ l/2 R cd h a c h b d + A" 1/2 ( 3 ) Rt a t h ] [(□ - ml) R ab + m\g ab R cd R cd \ = . (2.26) 

Recalling that t a t a = —A and noting that h ab R a b = ( - 3 ' ) R we see that Eq. ( 12.261) reduces to, 

[\ 1/2 R cd h a c h b d + \- 1/2 (3) Rt a t b ] (□ - mf) R ab = , (2.27) 

which is independent of the value of m 2 . By Lemma ( 111.21) . for mi ^ 0, this implies R ab = 
0, provided the space-time is static and asymptotically constant, whilst for mi = and 
the space-time being static and asymptotically flat implies R a b = 0. Indeed integrating 
Eq. ( 12.271) over the spatial slice exactly gives Eq. (I2.24j) . □ 

The final step that is required in order to allow us to consider Eq. ( 12. 9 p is to extend the 
above Lemmas to include a term of the form R ab R acb d. To do this note that by using the 
definition of the Riemann tensor, 

(V a Vfe — VfoV a ) R c d = R abc e R e d + RabdRce , (2.28) 

and the fact that the contracted Bianchi identity, for R = 0, gives, 

V c R b c = ^V b R = , (2.29) 

we can write (for R = 0), 

R cd Racbd = — V c V a -R b c + R ac R b ■ (2.30) 

Thus we see that this term is in fact closely related to OR ab and can be dealt with in a 
similar manner. 
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Lemma II. 4. For a static space-time the following relation holds: 
Vhd 3 x [(X^R^hlh) + A" 1/2 (3) Rt a t b ) R acbd R cd ] = 



x 



D a (-X 1/2 R ef D e Raf + A 1/2 (3) RD b R ab ) + A 1/2 (D a R ef ) (D e R af ) 



+A 1/2 (D a R ac ) (D b R bc ) - 2A 1/2 (D c R ca ) (D a (3) R) + A 1/2 \R a b R b c R c a - ( {3) R) 



(2.31) 



Proof. To show this we project Eq. (I2.30p onto the spatial surface, S, with h a e h K and use 
the identity 



h a e h b f V c V a R b c = h\h b f ( h cd - -t c t d ) V V /; 



c y a-n-bd 



1 

A 



1 ,( 3 ) „ 1 



D a D e R a f - — (D e X) (DfXy 3 ' R-— (D e X) {D d X) R fd 
+ 2T {DfX) ( ^ 3)jR) + ^ iD<1X) ( DeRfa) ' 



2A 



(2.32) 



where the second equality comes from the tedious but straight forward algebra coming from 
bringing the h b ^t d inside both covariant derivatives and repeatedly using Eq. ( 12. 4 p , Eq. ( 12. 5 j) , 
Eq. (12~T8D and Eq. ( 12191) . 

Substituting this expression into Eq. (I2.30p . multiplying by X 1 / 2 R e f and integrating over 
S, one finds, 



^/hd 3 x [X 1/2 R ef h a e h b R acbd R cd ] 



Vhd 3 x 



D a (-X 1/2 R ef D e R af ) + A 1/2 (D a R ef ) (D e R af ) + ^X~ 3/2 (3) RR ef (D e X) (D f X) 



+^X-' 3/2 R ef R fd (D e X) (D d X) - ^X- 1/2 R ef (DfX) (D {3) R) + X 1,2 R\R b c R c a 



(2.33) 



Similarly we project Eq. (I2.30p along t a t to find 



+a + b T>cd p 
tin tiadod 



—t a t h V(^aR b C + RacR b C t a t 

= i [(D a R a b ) (D b X) + R ab D a D b X +^ RD a D a X] - A [^R] 2 , (2.34) 

where the second equality follows from the staticity of the space-time and again, repeated 
use of Eq. ([H]), Eq. O , Eq.fl2HEJ) and Eq. (|2TT9|) . 



c+a+b 
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Multiplying Eq. ( 12.341) by A 1 ^ 2 ^ R and integrating over S, one finds, 
Vhd 3 x [X- 1 ' 2 ® Rt a t h R c dR acbd ] = 

Vhd 3 x D a Qa- 1 / 2 ® RR ab D b \ + ^X- 1 ' 2 ^ R^RD a \ 

+ 1 A -3/2 (3) ^ (DaA) {DbX) + 1 A -3/2 (3) ^(3 )i? 

_1 A -l/2 i? ab ^(3)^ (AA) _ A l/2 (3) £ ^(3)^ p« A) _ A l/2 J(3) fl ] = 



(2.35) 



Finally, adding Eq. (I2.33f) and Eq. (12. 35ft and using the contracted Bianchi identity given in 
Eq. (I2.25P demonstrates the desired result, Eq. (I2.3ip . □ 

III. GENERALISED LICHNEROWICZ THEOREM FOR 4 th ORDER GRAVITY 

Using the lemmas proved in the previous section, we will now extend the Lichnerowicz 
theorem to 4 th order gravity (for a^O and a > 3(3), provided the spatial scalar curvature 
satisfies certain bounds. 

Theorem III.l. Consider a static space-time, with a spatial slice that is topologically M?, 
in which the spatial curvature everywhere satisfies the following two conditions: 

m 2 - (3) R > , 
R a b R b a (m 2 + n) > , (3.1) 



pa pft pc 

R b R c R a 

R\R\ K } 



where 1Z is defined as, 



Let the space-time be asymptotically constant, unless both of the inequalities in Eq. \3. 1\) are 
saturated, in which case let the space-time be asymptotically flat. Then for such a space-time 
the equations R = and 

DR ab - m 2 R ab + ±g ab R cd R cd - 2R° d R acbd = , (3.3) 

imply R ab = 0. 
Proof. The equation, 

OR* - m 2 R ab + \g ab K d R cd - 2R° d R aM = , (3.4) 



implies, 



/ = {\ l l 2 R^h a e h b f + A" 1/2 (3) Rt a t b ) (nR ab - m 2 R ab + X -g ab R cA R cd - 2R cd R acbd ^j = . 

(3.5) 
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By Lemmas AO) and (UTil) and Eq. f[2T2D and Eq. ( I2~25j) of Lemma fllL2]) the integral of 
Eq. (13. 5p over the constant time slice, S is, 

J y/hd 3 x (I) = J \fhd?x^D a [X 1/2 (^RD^R + R cd D a R cd - R ef D e R af RD 6 ^)] 
-A 1/2 (D a R cd ) (D a R cd ) - 2A 1/2 (D a ^ 6c ) (D b R ac ) - A 1/2 [D a < 3 > i? - 2D c R a c ] [D® _ 2 J D 6 i2 , 



-A 1/2 \R cd R cd (m 2 + ft) + (^i?) 2 (m 2 -® R) 



. 



(3. 



Asymptotic constancy ensures that the boundary terms vanish, whilst the final two terms 
are negative definite if the spatial curvature satisfies the properties, 



m 2 - (3) R > , 
R\R h a (m 2 + n) > . 

Finally, as shown in Appendix AVID , the combination 

(D a R cd ) (D a R cd ) + 2 (D a R bc ) (D b R ac ) > 



(3.7) 



(3.8) 



for all R a b. Thus each term in the integrand of Eq. (13. 6p is negative definite and hence 
required to vanish. Thus D c R a fy = and if either of the inequalities in the expressions 
given in Eq. (13. 7p fail to be saturated in any open region, Eq. (I3.6P implies R a b = 0. If the 
expressions in Eq. (13. 7p are saturated everywhere, then asymptotic flatness is required to 
imply R ab = everywhere. □ 

The equations of motion of 4 th order gravity (with a — 3/3 > 0), imply R = (by 
Theorem (jTLT]) ) and (from Eq. RZ^ ) 



ha ( DR ab - m z R ab + -g ab R cd R cd - 2R cd R acbd ] = 



(3.9) 



with m = 7/t h a . Thus, provided the spatial curvature everywhere obeys the bounds, 



Mr < -J- 

7 



n > 



han 2 



(3.10) 



there are no, non-trivial, static, asymptotically constant vacuum solutions to 4 th order grav- 
ity (with a — 3(3 > and q/0). 



IV. GENERALISED ISRAEL THEOREM FOR 4 th ORDER GRAVITY 

The result proved in Section demonstrates that there are no static, vacuum solutions 
to 4 th order gravity for space-times which have a spatial section that is topologically M 3 , but 
one may ask whether a similar result holds for space-times in which the spatial slice contains 
a boundary. In particular, in order to consider space-times with a (static) black-hole, we 
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need to allow for a null, interior boundary to our space-time i.e. the equal time hyper-surfaces 
become null. At first sight the presence of boundary terms in Eq. (13. 6p would appear to be a 
significant difficulty for such a space-time, however on the horizon, the normal to the spatial 
hyper-surface, n a becomes null i.e. n a n a = A = and hence the boundary terms vanish. 
Thus, if the space-time is asymptotically constant (which ensures the boundary terms at 
infinity also vanish) and all interior boundaries are null, the result of Section (IIIII) will still 
apply. 

Theorem IV. 1. Consider a static space-time, in which the spatial curvature everywhere 
satisfies the following two conditions: 
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hOLK 2 



n > --^ , (4.i) 

Let the space-time be asymptotically constant, unless both of the inequalities in Eq. \4-l\j 
are saturated, in which case let the space-time be asymptotically flat. Let the spatial slice be 
bounded by a null surface that is topologically S 2 . Then the only solution to 4 th order gravity 
(with q^O and a — 3/3 > 0), in the region exterior to the null surface, is R a b = 0. 

Proof. As in Theorem lIII.lt the equations of motion for 4 th order gravity imply, (Eq. (13. 6p ) 

J Vhd 3 x (J) = J Vhd 3 xl^D a [A 1/2 (®RD®R + R cd D a R cd - R ef D e R af -® RD b R ba )] 
-A 1/2 (D a R cd ) (D a R cd ) - 2A 1/2 (D a R bc ) (D b R ac ) - A 1/2 [D a (3) R - 2D c R a c ] [D® - 2D b R ab ] 
-A 1 / 2 [R cd R cd (^a- 1 + TZ) + (^R) 2 ^oT 1 R)] | = . (4.2) 

At infinity the boundary terms vanish because the spatial section is asymptotically con- 
stant, whilst on the interior boundary they vanish because A = 0. Thus, just as in Theo- 
rem (IIII.ll) . provided the 3-dimensional curvature satisfies the conditions 
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- (3) R > 



han 2 

each term in the integrand of Eq. (14. 2 j) is negative definite and hence vanishes. If either of 
these inequalities is not saturated everywhere, asymptotic constancy is sufficient to imply 
R a b = 0, otherwise asymptotic flatness is required. □ 

Since 4 th order gravity implies R ab = (for asymptotically constant and static space- 
times), provided 
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even in the presence of null boundaries, it follows that the only spherically symmetric solution 
is the (exterior of the) Schwarzschild metric. 

It is important to note here that this result relies on the existence of a null boundary to 
the equal time hyper-surfaces and the existence of a vacuum (T ab = 0) everywhere between 
this horizon and infinity. This is true for a (static, vacuum) space-time containing an event 
horizon, however the result cannot be applied to the exterior of a general spherical object that 
does not contain such a horizon. In particular this result to does not imply that the metric 
outside a spherically symmetric distribution of matter will be Schwarzschild, unless that 
matter is entirely contained within the event horizon. This agrees with the (perturbative) 
results of (for example) (6|. 



V. NON- VACUUM SOLUTIONS 



In Sections (IHip and (lIVp we showed that for T ab = 0, the static, asymptotically constant 
solutions of GR are also solutions to 4 th order gravity (with a^O and a — 3(3 > 0), at least 
in the case of space-times that have a spatial topology of M 3 (Theorem (1III.1I) ) or space-times 
that have internal null boundaries (Theorem ( llV.lj) ). and whose spatial curvature everywhere 
satisfies Eq. (14. 3p . A natural question that one may ask is whether this connection between 
the solutions of GR and 4 th order gravity extends to the non-vacuum case i.e. T ab ^ 0. We 
will now demonstrate that it does not. In fact below we show that for, static, asymptotically 
flat space-times containing a barotropic fluid with equation of state to > — ~ —0.6933, 
only the vacuum solutions of the two theories coincide, with all other solutions being distinct. 

Theorem V.l. Consider an asymptotically flat solution to GR, GR g a b, with matter satisfy- 
ing, 

Vhd 3 xT a b T b c T c a > . (5.1) 



Then GR g ab is a solution to 4 th order gravity (with a^0 and a ^ 3(3) iff R ab = 0. 
Proof. Consider a solution to the Einstein's equations: 

7«T 2 (^R ab - ^gabR^j = l,T ab , (5.2) 

that is static and asymptotically flat. Substituting this into the equations of motion for 4 th 
order gravity, Eq. (12. ip we find, 

{a - 2(3) R. a . b - aUR ab 4/3) g ab DR + 2aR cd R acbd 

-2(3RR ab - ^g ab (aR cd R cd - (3R 2 ) = . (5.3) 
The trace of Eq. (15.31) is, 

(a — 3(3) OR = . (5.4) 
Then by Theorem ( III. Ill , for a — 3(3 ^ 0, this implies R = 0. Thus Eq. ( 15. 3p becomes, 



— a 



DR ab + -g a bR cd Rcd — 2R cd R acbd 



, (5.5) 
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which by Theorem (IHI.lj) implies (see Eq. ( 13.6 



J Vhd 3 x (I) = J Vhd 3 x!^D a [A 1/2 (®RD®R + R cd D a R cd - R ef D e R af -® RD b R ba )] 
-A 1/2 (D a R cd ) (D a R cd ) - 2A 1/2 (D a R bc ) (D b R ac ) - A 1/2 [D a < 3 > R - 2D c R a c ] [D® - 2D b R 



-A 1 / 2 \R a b R b c R c a - (Wr) : 



. 



(5. 



Thus, if 



Vhd 3 x \R\R b c R c a - ( (3) i?) : 



> 



(5.7) 



then each term in the integrand is negative definite and hence vanishes. 
However it is easy to show that, 



pa pb pc 
K b K c K a 



fa j 



g ad g be g cf RbdRceR 

R\R b c R c a -{^R) S , 



(5.8) 



h ab — jt a t b , using the Codacci 



where the second equality follows from decomposing g ab 

and Gauss equations (Eq. (I2.19p and Eq. (j2.18|) ) and noting that R = 0. Thus each term in 
the integrand of Eq. (I5.6P is negative definite, and hence vanishes, provided 



Vhd 3 R a h R b M c „ > . 



(5.9) 



Using Eq. (15. 2\i and the fact that R = for these solutions, this condition can be written as 

Vhd 3 xT a b T b c T c a > . (5.10) 



Thus, for matter satisfying Eq. (I5.10p . a static, asymptotically flat solution to GR is also 
a solution to 4 th order gravity iff R a b = 0, which, by Theorems (IIII.ll) and ( 1IV.1I) corresponds 
to the (static, asymptotically flat) vacuum solutions in both GR and 4 th order gravity. Thus, 
other than the vacuum, there are no static, asymptotically flat solutions to GR, with matter 
satisfying Eq. (I5.10p . that are also solutions to 4 th order gravity. □ 

If we consider a perfect fluid the energy momentum tensor can be written as T a b = 
pu a Ub + P (S a b + u a Ub), with u a a unit time-like vector tangent to the observer's world-line, 
P the fluid's pressure and p its energy density. Eq. (I5.10p then becomes, 



J Vhd 3 x (p 3 + 3P 3 ) > , 



(5.11) 



Assuming the strong energy condition holds, p > 0, then this condition is trivial for all 
matter with positive pressure. If we consider a barotropic fluid with equation of state 
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parameter u> i.e. P = up then Eq. (15. 10p reduces to, 



J Vhd 3 x (1 + 3u 3 ) > 



(5.12) 



This is, in particular, satisfied for tu > — 3 -1 / 3 ~ —0.6933 everywhere, and hence for dust 
(u> = 0) and radiation (u = 1/3) fluids. 

Thus for a static, asymptotically flat, space-time containing a single, barotropic fluid for 
which u > — 3 -1 / 3 everywhere, no solutions to GR remain (static and asymptotically flat) 
solutions to 4 th order gravity. 

One can also consider, for example, the energy-momentum tensor of a Maxwell field, 



T\ = — ( F a r F c h ' 
b 4vr V c b 



- A 9\F d e F e d 



(5.13) 



with F a b the electromagnetic field tensor. For such a field, one can directly check that 

1 



rpa rpb rpC 



64tt 3 



I (TrF 2 ) 3 - ~ (TrF 4 ) (TrF 2 ) - TrF e 



> 



(5.14) 



where the matrix F is given by F a b and hence is anti-symmetric and the final inequality 
follows from direct computation, for a general anti-symmetric matrix. 



VI. CONCLUSIONS 

In the Sections (1IIII) and (llVj) we demonstrated that there are no, non-trivial, static, 
asymptotically constant solutions to 4 th order gravity, provided the 3-dimensional curvature 
satisfies the following two conditions: 

®R < and K > -—^ , (6.1) 

to fir naK z 

everywhere outside a null horizon. The meaning of these conditions is the following. We 
require that the 3-dimensional scalar curvature, measured by both and 7Z, to be ev- 
erywhere smaller in magnitude than the scale set by the ratio between r yn~ 2 and ha. This 
is exactly the scale at which corrections to General Relativity will become significant. If 
these corrections are motivated by quantum corrections to general relativity then this will 
be the quantum gravity scale. Thus our results can interpreted as follows: there are no, 
non-trivial, static, asymptotically constant, vacuum solutions to 4 th order gravity, so long 
as the 3-dimensional scalar curvature is everywhere less than the quantum gravity scale. In 
particular for black holes, this implies that the Schwarzschild solution is the unique, spheri- 
cally symmetric, static, asymptotically constant vacuum solution to 4 th order gravity, unless 
the 3-dimensional scalar curvature exceeds the quantum gravity scale outside the horizon, a 
condition that is easily achieved for macroscopic black-holes. 

This also demonstrates the fact that the 'no-hair' theorem applies to 4 th order gravity 
provided that the spatial curvature (outside the horizon) is everywhere less than the quantum 
gravity scale. Note however that the converse does not follow i.e. if the spatial curvature 
outside the horizon is greater than the quantum gravity scale, Theorem (lIV.ip does not 
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imply that there are additional (static, asymptotically flat) solutions. 

Finally, in Section (jVj), we demonstrated that for matter satisfying the rather mild con- 
dition, 

J Vhd 3 xT a b T b c T c a > , (6.2) 

except for the vacuum, none of the static, asymptotically flat solutions to GR are solutions to 
4 th order gravity. In particular this includes space-times containing single, barotropic fluids 
with equations of state everywhere satisfying, u > — 3 -1 ^ 3 , such as relativistic and non- 
relativistic Baryons and Fermions. A specific application of this is that the GR solutions 
for (static, asymptotically flat) stars will fail to be solutions in 4 th order gravity (for the 
same matter distribution). This agrees with the results of [6], which gives the solutions 
to 4 th order gravity for a particular, extended, spherically symmetric, system (calculated 
perturbatively and close to the origin) which is shown not to be Schwarzschild. 
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Appendix 

In this appendix we prove various useful (if rather laborious) identities. First consider 
t a 'V a t b 'VbRcd, for a static space-time. Using the staticity condition, 

£>tRab = o , 

t a V a R cd = -2R a[d V c) t a , (6.3) 

we can write, 

t a V a t b V b R cd = 2 [(V b t a ) R a(d (V c) t b ) + R ab (V {d t a ) (V c) t 6 )] 

-R bd t a W a W c t b - R ch t a V a V d t h . (6.4) 

Using t c V c A = and hence that V a (t c V c A) = 0, we find 

t c V c V a A = - (V a t c ) (V C A) . (6.5) 

Thus, 

t c V c V a t d = ±- (t c t d V c V a X - t a t c v c v d x) 

= ^ H d (V a t c ) (V C A) + t a {V d t c ) (V C A)) = , (6.6) 
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where the final equality follows from expanding V a t c = j- (t[ c V a ]A). Using this in Eq. (16 .4p 
gives, 

t a V a t b V b R cd = 2 [(V b t a ) (R <d V c) t b ) + R ab (V c t b ) (V d t a )} , (6.7) 

which is Eq. f[2TT0]) . 

By using V a t c = v (£[ c V a ]A) and recalling that t a V a A = 0, we find, 

(V b t c ) (V a t d ) = ^ (t c V,A - t d V c \) {t d V a \ - t a V d \) , 

= ^ 2 t c t a (V d A) (V d A) + 1. ( V C A) (V«A) , (6.8) 

which is just Eq. (12. lip . 

In order to show Eq. (I2.14p note that, 

jt c t d V b R cd = V b {^t c t d R cd ^j - R cd V . (6.9) 

Now we use the fact that for a static space time, the extrinsic curvature of the constant 
time hyper- surfaces vanishes (this can be checked by direct computation). Thus the Gauss 
equation (Eq. (I2.19P ) gives, 

\t a t b R ab = (3) R- R , (6.10) 
A 

whilst the Codacci equation, (Eq. ( 12. 18[) ) is simply t a h b c R ab = 0. A consequence of the 
Codacci equation and t c V c A = is, 

t a (V b X)R ab = , (6.11) 
Using this and the Gauss equation in Eq. (I6.9P gives, 

jt c t d V b R cd = V b { (3) R-R) . (6.12) 
With this we directly find Eq. fl2TT4"|) : 

^h ab (V a A) t c t d V b R cd = ^ (V a A) (V a ( {3) R - R)) . (6.13) 
Similarly, by noting, 

1 



jt c t d h ab V a h e b V e R cd = h%V a \^t c t d h e b V e R cd ) - h ab 
and using Eq. ( 16.121) we find, 



V„ ( -t c t d 



V b R cd , (6.14) 



jt c t d h ab V a h e b V e R cd = D a D a (®R - R) , (6.15) 

which is Eq. (12T7p . 

A key equation of Lemma f ill. 2h is Eq. (12 . 1 6j) . which we derive here. By decomposing 
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9ab = Kb ~ jt a t b we have, 

DR ab - m 2 R ab = h ab V a V b R cd - jt a t b V a V b R cd - m 2 R cd = . (6.16) 

Splitting the projection h ab = h ac h b c and bringing one of the 3-metrics, inside the covariant 
derivative the first term becomes 

h ab V a V b R cd = h ab V a h e b V e R cd , (6.17) 

where we used the fact that because V a Qbc = 0, we have 

V a h bc = V a (jhtc) , (6.18) 

and h ab, Vatb = by symmetry. For the second term in Eq. (16.161) . we again use the fact that 
VcSu = to find, 

V c /i a b = V c ( —t a t b 

t a V c t a = -t a V a t c = -^V C A , (6.19) 

where the second line follows by expanding the right hand side, taking the trace over (a, b) 
and noting that V( £&) = 0. Using Eqs. (16. 7\i and (16.171) in Eq. (I6.16P gives, 

h ab V a h\V e R cd + ^h ab (V B A) (V b R cd ) 
+^ [t c t a (V 6 A) (V 6 A) - A (V C A) (V a A)] R ad 
-2^2 (V d A) - t d (V a A)] [t b (V C A) - t c (V d A)] R ba 

+ lp ( V * A ) ( V ' A ) " A ( V ° A ) ( V ^ A )] R ca ~ ™ 2 Rcd = , (6.20) 
which is Eq. flgjgfl . 

Finally, we need to show that the terms given in Eq. ( 13. 8 p are positive i.e. that 

(D a R bc ) (D a R bc ) + 2 (D a R bc ) (D b R ac ) > . (6.21) 

To do this, consider an arbitary tensor T abc that is symmetric in the final two indices i.e. 
T a bc = T acb . Now we use a tetrad basis, e a , that is othonormal i.e. e % a e°- = Sp to decompose 
the (positive definite) metric as, 

Kb = j a 4 S U ■ ( 6 - 22 ) 



With this one can write 



J-abc = II n rl J-abcJ-de} = J-ijkJ-mnpO J y , 

T abc T bac = h ae h bd h cf T abc T def = T ijk T mnp 5 m S jm 5 kp . (6.23) 
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We now form the combination, 

T abc T abc + 2T abc T bac , (6.24) 

and expand out the Einstein summations. Since T abc is a three index tensor, the three 
possibilities we need to consider are: the indices are the same, one is distinct and all are 
distinct. The first possibility clearly involves only positive terms. In the orthonormal tetrad 
basis given by Eq. (I6.22p these are 

T abc T abc + 2T abc T bac = 3 (T m T m + T 222 T 222 + T 333 T 333 ) + ... , (6.25) 

where the dots remind us that there are other terms we need to account for. The terms in 
Eq. f)6.24p that arise from one index being distinct, with the other two being the same, can, 
by using Eq. (I6.22p . be grouped to give, 

T abc T abc + 2T abc T bac = (2% + Tijj) (2T Uj + T w ) + . . . , (6.26) 

where there is no implicit summation over repeated indices and again the dots indicate that 
there we are considering only certain terms. Finally we need to account for the terms in 
Eq. (16.24p for which all the indices are distinct, which are given by, 

T abc T abc + 2T abc T bac = 2 (Tijk + Tj ik + T ki j) (T^ + Tj ik + T ki j) + . . . , (6.27) 

where again there is no implicit summation over repeated indices and the dots indicate that 
there are additional terms. In Eq. (16.261) and Eq. (I6.27p . crucial use has been made of the 
fact that Tij k = T ik j. Finally, putting these together we find that, in this orthonormal tetrad 
basis, 

T T abc -\- 2T T bac 3 ^ ^ TnjTiii -\- ^ ^ {2Taj -\- Tijj) {2Taj -\- Tijj) 

i i^j 

+ 2 {Tij k + Tji k + T^j) {Tijk + Tjik + T^j) , (6.28) 

where once again there is no implicit summation over repeated indices. Clearly each term 
in Eq. (I6.28j) is positive. It is worth noting here that this result hold in any number of 
dimensions, provided the metric, h ab , is positive definite. 

One may be concerned by the fact that if T abc were also anti-symmetric in the first 
two indices, the combination given on the left-hand-side of Eq. (16.281) . would appear to be 
negative. However there are no three index tensors that are anti-symmetric in the first pair 
of indices and symmetric in the second pair, as can easily be checked, 

T a bc T bac T bca T cba T cab T acb T abc . (6.29) 

Eq. (I6.28P holds for any three index tensor that is symmetric in its final pair of indices 
and in particular it holds for T abc = D a R bc , thus, 

(D a R bc ) (D a R bc ) + 2 (D a R bc ) (D b R ac ) > , (6.30) 
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as required. 
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